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O^. Abstract 
Vh . 

Evolution of the density parameter in the anisotropic DGP braneworld model is 
studied. The role of shear and cross-over scale in the evolution of fip is examined 
^ ■ for both the branches of solution in the DGP model. The evolution is modified 

0^ , significantly compared to the FRW model and further it does not depend on the 

I value of 7 alone. Behaviour of the cosmological density parameter Qp is unaltered 

•/^ I in the late universe. The study of decceleration parameter shows that the entry 

CN ' of the universe into self accelerating phase is determined by the value of shear. 

We also obtain an estimate of the shear parameter ~ 1.68 x lO^^'^, which is in 
agreement with the constraints obtained in the literature using data. 

^ : 1 Introduction 

Results from the recent observations of the type-I supernovae and WMAP 
suggest that we are hving in more or less flat type universe [H [21 El S] • These 
implies that the cosmological density parameter Q must be equal to unity 
or very much near to it. The inflationary scenario also supports the value 
of the density parameter which is very much close to unity. For the esti- 
mation of the cosmological density parameter to be unity it is assumed that 
apart from the dark energy and baryonic matter the universe also contains 



*Eniail: rizwan@iucaa.ernet.in 
''^Email: pkssp@uohyd.ernet.in 



1 



dark matter. Therefore the contents of the universe are of paramount rele- 
vance in understanding the geometry and evolution of the universe. Hence 
determination of the density parameter is very important because it can put 
constraints on cosmological model building and that are subjected to obser- 
vations. In general the density parameter can be considered as dynamical 
quantity rather than a constant. For the evolutionary study of the density 
parameter one has to take into account the role of the effective equation of 
the state of the epoch of the universe under consideration. Most of the recent 
observations estimate the value of total density parameter Qq to be close to 
unity, assuming the homogeneity and isotropy of the spacetime. However, 
one can also consider the universe that deviates from such idealized situation, 
to accommodate an anisotropic spacetime. Historically, anisotropic models 
of the universe have been studied to avoid assumptions of specific initial 
conditions in the Friedmann-Robertson- Walker (FRW) models and also it is 
believed that dynamics of early universe is profoundly influenced in presence 
of anisotropy just below the Planck or string scale O [6]. Bianchi models are 
also relevant in the context of Belinskii-Khalatnikov-Lifshitz (BKL) conjec- 
ture [7] which states that in the Einstein equations 'terms containing time 
derivatives will dominate over those containing spatial derivatives' as one 
approaches a space-like singularity, such a model is described by homoge- 
neous and anisotropic cosmologies. Evolution of the cosmological density 
parameter in presence of anisotropy is studied and showed that the initial 
conditions are robust with respect to shear and viscosity [S]. In the present 
work our aim is to examine dynamical behaviour of the density parameter 
in the anisotropic Dvali-Gabadadze-Porrati (DGP) braneworld cosmology, 
which has the potential to explain the current accelerating phase of the uni- 
verse. The accelerated expansion of the universe is again suggested by the 
recent observations such as type-I supernovae and WMAP[T1 El El II]- The 
DGP braneworld cosmological model has also been studied from observa- 
tional point of view and these have put constraints on model parameters 

The DGP braneworld cosmological model is basically inspired by higher 
dimensional theories that have been studied in great detail in recent times. 
In these models our observed four dimensional (4D) universe is a three di- 
mensional hypersurface known as brane embedded in a higher dimensional 
spacetime called bulk (for a review of brane cosmology see [H] ) . The appli- 
cation of such scenario to the homogeneous and isotropic FRW brane leads 
to modification of the Friedmann equations with a quadratic correction to 
energy density at higher energies. Many issues in cosmology like inflation, 
dark energy, cosmological constant were investigated in the braneworld cos- 
mological scenario and interesting results were obtained. Later the studies 
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were extended to anisotropic brane, such as cosmological solutions to the 
Bianchi-I and V braneworld models were obtained and it is shown that for 
matter obeying barotropic equation of state, there could be an intermediate 
phase of high anisotropy which in turn leads to an isotropic phase in the later 
universe [12]. Another interesting result is obtained in [13], with the scalar 
field as source, a large initial anisotropy induces significant damping in the 
dynamics of the scalar field resulting in greater inflation. The shear dynam- 
ics on Bianchi I brane cosmological model with a perfect fluid is considered 
and shown that shear parameter has a maximum during a transition period 
from non-standard to standard cosmology [H]. 

Motivated by these developments, we consider evolution of the density pa- 
rameter in an anisotropic DGP braneworld model and also examine the late 
time acceleration in this model. Anisotropic brane in the DGP model is also 
considered in [15l |16] and solutions are obtained without affecting the basic 
features of the model. The DGP cosmological model possesses two classes of 
solutions; one which is close to the standard FRW cosmology and the other 
is either a fully five dimensional regime or a self-inflationary solution which 
produces accelerated expansion. The self-accelerated branch is known to con- 
tain ghost and tachyonic like excitations [T71 [181 US EHl [21] , but nevertheless 
its ability to explain the late time acceleration makes it an interesting case 
for study and has generated enormous amount of interest in the recent years. 
We mainly focus on the effects of shear and cross over scale on the dynam- 
ics of dimensionless density parameter and the deceleration parameter. The 
evolution equation for Qp is derived in the self-accelerated branch and for the 
fully five dimensional branch and the corresponding deceleration parameter 
are also obtained. 



2 DGP dynamics on Bianchi-I brane 

In the DGP model our universe is a 3-brane embedded in five- dimensional 
(5D) Minkowski bulk and there is an induced four- dimensional Ricci scalar 
on the brane, due to radiative correction to the graviton propagator on the 
brane [2^. In this model gravity is modified at large length scales and there is 
a crossover scale between the four- dimensional and five-dimensional gravity 
given by 



where and //^ are the 4D and 5D gravitational constants respectively. 
Below the crossover scale Tc, the potential has the usual Newtonian form and 
above which the gravity becomes five-dimensional (A similar scenario with 
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a scalar curvature term in the brane action is considered in [251 with a 
brane embeded in AdS bulk. Here, it is shown that it is possible to generate 
a four-dimensional theory of gravity depending upon AdS lengths). We start 
with a generalized DGP model in which both bulk cosmological constant A 
and brane tension a are non-zero. The effective Einstein equation on the 
brane can be written as 

+ ^rf,^ + -\F^u + k^U^u + -^L^^ - E^^, (1) 

n^.. = -^r^pr(^ + ^rr^u + ^q^,uTa^T'''^ - ^Q^ut'^. (2) 
Pnu = —-^G^pGI + —GG^u + -QfiuGa/sG"'^ — -^q^iuG'^, (3) 



L,. = \{G,,r^ + T,pG^,) - ^{tG,. + Gr,,) - \q,.{G^pr^^ - \gt), (4) 

where t^^ is the energy-momentum tensor of the matter fields on the brane, 
a is the brane tension and T^i, is the bulk energy-momentum tensor. It can 
be seen that apart from the quadratic matter field corrections Tl^v to energy- 
momentum, there are corrections coming from the induced curvature term 
through F^ij and L^^. E^^, are the non-local corrections, coming from extra 
dimensions which correspond to the projection of the bulk Weyl tensor on 
the brane. If we define the four velocity comoving with matter as u^, the 
non-local term takes the following form, 

E^u = -^rl[U {u^u^ + ^h^^) + P^^ + Q^Uy + QyU^]. (5) 



Where h^^ = g^^ + u^u^ and 



U = -l-Ep,uV, (6) 



c 

is the effective non-local energy density on the brane 



c 



- Ih'^^hJ E^p (7) 
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is the effective non-local anisotropic stress and, 

= ^K'^E^pu^ , (8) 

c 

represents the effective non-local energy ffux on the brane. The brane en- 
ergy momentum satisfies the conservation equations V^r^j, = 0, the detailed 
description of conservation and dynamical equations is given in |25j . 

We are interested in anisotropic cosmological models on the brane, hence 
we consider the Bianchi I brane model. We start with the 5D metric of the 
following form 

ds'^ = dy'^ + qf.jydx^dx", (9) 

where the brane is located at y = 0. The Bianchi metric on the brane is 
given by, 

rfs\=o = q^u{y = ^)dx^dx'' = -dt^ + a'f{t)dx'^ . (10) 
The conservation equation takes the following form, 

p + e(p + p) = o (11) 

D'P^, = (12) 

U + IQU + (T^'^^P^, = (13) 

where dot denotes the u'^V^ and G represents the volume expansion rate 
and a^'^ is the shear rate. There is no evolution equation for P^,y, since it 
is bulk degree of freedom and cannot be predicted from the brane [25]. The 
Hubble parameters for the Bianchi I metric is given hj Hi = ^ and one can 
define the mean expansion factor as = (010203)^/^, thus 

e^3H = 3^^Y.H^■ (14) 

The modified Raychaudhuri equation for the Bianchi-I metric on the 
brane is obtained as [15], 

2 



l„o „.. -k 



2 , p. 2 



K\ o \^ 1 



^[H' + — -fi'p +-p(p + 2p) 



2p2 V V 52 



5) 
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and Gauss-Codacci eqns are, 



^t,v + 00-^.. = ^KP^u (16) 

+ a^^a,. = |t (^3 [h^ + ^) - ^?p) ' + ArlU . (17) 

The system of equations are not closed on the brane due to the presence 
of anisotropic stress term P^,^ and complete bulk solutions are needed to 
determine brane dynamics. This problem is avoided(as in earlier papers 
[T3| [T5| |25] ) by choosing a special case U = 0, this vanishing of non-local 
energy density leads via f lT3|) to condition a^uP^^ = 0. This consistency 




condition implies a condition on evolution of P^u as cr^'^P^i, = Ar'^P'^'^P^,^, 
following from Eq. f lTB]) . As there is no evolution equation for P^,^ on the 
brane, this is consistent on the brane (for complete discussion see [121 125]). 
This allows us to find the dynamics of shear scalar, equation ( 1T6|) can be 
contracted with shear and integrated to get, 

a^V^, = 2a2 = — , S = 0. (18) 

The modified Friedmann type equation in the DGP cosmology for the 
Bianchi-I metric flTUl) is obtained using eqns f[T7l) and flTSl) as 

(^^ + f ) - 1^ - (^^ + ^) - ^ = (19) 

where e = ±1 corresponds to two possible embedding of the brane in the 
bulk space-time and gives rise to the two branches of solutions in the late 
universe. We can rewrite f[T^ as 



(20) 

(21) 
(22) 

(23) 
(24) 
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The Friedmann type equation fl2U]) takes the form 



1 = (^y^+yoT+^^j +^^E + ^^x (25) 

which imphes the variables must take values in the interval [0,1] and further 
it is seen that the normalisation condition is modified from the usual FRW 
case. For the K = isotropic DGP model, people have estimated values of 
fire compared with observation and which are in good agreement 

with theory |26]. For the Bianchi I case we have to take into account Qy. to 
get correct estimate of values of ^2,.^ and flp, as can be seen from fl25l) . except 
in the limit 5* — t- oo when effects of shear are negligible. 




0.2 0.4 0.6 0.8 1.0 ' 

Figure 1. Evolution of Qp for e = +1, for different values of fipQ, with S = 0.001, 

A- = 0,7 = 4/3 and = 0.1. 



3 Evolution of cosmological density param- 
eter 

Next we study the evolution of density parameter Qp and see its dependence 
on shear a and cross over scale Vc- For this we consider the two branches of 
solutions (e = ±1) of the DGP scenario separately. The branch e = +1 cor- 
responds to the self-accelerating solution in the late universe, without a cos- 
mological constant. Expanding equation (120|) under the condition /i^pr^ ^ 1 
and considering e = +1 case up to first order |T51 127], we obtain the Fried- 
mann equation of the form, 

9 ^ 1 2 , E2 



Using the definition fl2T|) the above equation can be written as, 

^'(2^p - 1) = - §^ + f ' (27) 
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it can be noted from the above equation that K = no longer imphes that 
fip = 1, as in the standard FRW case, we see the extra dependence on 
which is a typical feature of the DGP model. Also the effects of shear S 
cannot be neglected and value of flp is constrained by these parameters. We 
also derive the following quantity and will use it later 

H = -H'- + P) + l^^+ ^/iV - 2|^- (28) 




0.2 0.4 0.6 0.8 1.0 



Figure 2. Evolution of Vtp for e = +1 in the isotropic case ( S = 0) and for different 
values of S, with Op^ = 0.04, = 0, 7 = 4/3 and = 0.1. 

In order to study the evolution of fip, we differentiate eqn fl2T]) and use 
eqns ([26]) and ([28]) to get. 



dVLp 
dr 



HQ, 



(37-2)(2fip-l) 



1 2 



1 4S^ 



(29) 



Where we have assumed that the total energy density and pressure are 
related by linear barotropic equation of state, p = (7 — l)p. Next, to obtain 
equation for the Qp{S) phase plane analysis we divide eqn f l29l) by ^ = HS 
and substituting H from fl27|) we get. 



dQr 



Vtr. 



(37-2)(2^], 



(2^]p 



dS s 

The above equation can be integrated to get. 



_2_ I 4S^ 

_ J S£ I _^ 

r2 g6 ~r c2 



S2 



(30) 



1 



Qp, SI - 3KSy, + SV2 



So J 



37-6^ 



(31) 



here flpg and Sq denote the values of density parameter and scale factor at 
the present time. It is seen from (1311) that evolution equation depends on 
the shear and its behaviour is significantly effected in the early universe, 
depending on the value of S. And the value of f2p as S* — )■ is around 0.5 for 
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all values of 7 in the range (0,2), whereas in the FRW it depends on the 7 
and equation of state. We analyse this further by giving phase plane plots of 
flp, Figure. 1 shows behaviour of flp as a function of S for different values of 
Qpg in the self accelerated branch of the solution. The value of Qp for S very 
close to zero is same for different values of Qp^ and starts deviating for large 
values of S. Figure 2 shows the comparison of fip in the isotropic case along 
with different values of S in the anisotropic case. It is observed that for the 
isotropic case Qp starts with a high value, whereas for anisotropic case i.e. 
E = 0.1, 1 the plot is peaked at intermediate value of S and starts decreasing 
for large values of S. The dependence of density parameter on 7 is shown in 
the figure 3, the asymptotic value of fip as 5* — )■ is around 0.5 irrespective 
of the value of 7. And the density parmater decays faster for small values 
of 7. Finally the effects of cross over scale on the evolution of cosmological 
density parameter is shown in Figure 4. We have shown that changing the 
value of crossover scale i.e. Tc = 0.1, 1, 2 does not vary Qp much, whereas if 
we simultaneously change the value shear parameter S to a higher value, the 
behaviour is significantly modified for the same values of r^. 

-7=1 
y = 4/3 
— 7 = 2 



Figure 3. Evolution of Qp for e = +1 with different values of 7 and = 0.04, 
K = 0,1: = 0.001, rc = 0.1. 

We can define the following dimensionless deceleration parameter q = 
— ^ or q = to study the accelerated phase of the universe, which can 

be written using (128|) as, 

,= (3.-2H-(2n,-l)-^,-^,^^. (32) 

Acceleration takes place when g < 0, which is determined in the standard 
FRW case by the value of 7 and 7 = 2/3 forms the critical value. It is clear 
from equation fl32|) in this case entry into the accelerating phase is determined 

-'^For plotting the isotropic case we have used the expression Hp = 
1 

Opo S^-if2,.^ \, So ) 
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by values of the shear and cross-over scale. Figures 5 shows the behaviour 
of deceleration parameter as a function of scale factor, Fig. 5 (a) shows the 
dependance of q on the shear and it is seen that higher value of the shear 
imply that the universe enters an accelerated phase at later time. In figure 
5(b) the decceleration parameter is plotted for different values of rc, it is 
seen the irrespective of the value of cross over scale, the universe enters a self 
accelerated phase at the same time. And in both the cases (5(a) and 5(b)) 
it is noted that in the limit S* — )■ 0, g approches a value 2. 




i ' (b) 




Figure 4. 



Evolution of ilp for e 



PO 



0.04, K 



+ 1 for different values of Vc and 7 = 4/3, 
0. In the left panel we have taken the value of S = 0.001 and in 



the right panel S = 0.01, to show that a higher value of shear will lead to large 



variation in VLp with Tc 



Next, we consider the e = — 1 case (which is also known as the fully 
5D case [2H]), again expanding (I2U]) with condition /i^pr^ ^ 1 we get the 
Friedmann type equation as. 



K 

52 



Ac 2 ^2 

9 ^ 56 



(33) 



Since it represents the fully 5D case, the Freidmann equation has only quadratic 
contribution from the energy density. We define new dimensionless parame- 
ter as 



this allows us to rewrite the equation (15^ in the form 



(34) 



H^{np2 - 1) 



52 ~ 



and 



(35) 
(36) 
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Next, we discuss evolution equation for the dimensionless density param- 
eter, for this we follow the method used in e = +1 case. Differentiating 
with respect to time and using equations (135|) and (l36i) we obtain, 



dVLp2 
dr 



H{np2 -i)np2 



2(37 - 1) + 



(37) 



Integrating equation (!37j) gives 



l-VL 



KS^ 




(3^ 



where Vtp2 is the value at the present time. Let us analyse the effects of shear 
on the behaviour of VLp2. Fig.6(a) shows evolution of VLp2 for the the e = — 1 
case, it is seen that asymptotic value of VLp2 for 5 — t- is around 1 for all 
values of VLp2. But, if we take different values of 7 the behaviour changes 
drastically, figure 6(b) shows VLp2 for different values of 7 for K = {]. It is 
observed that for 7 = 4/3, 2 in the limit S* — )■ the value of VLp2 is around 1 
and decays from there. On the other hand for 7 = 2/3 in the limit 5* ^ 
the value of VLp2 — and increases for large values of S. 

Now the deceleration parameter can be obtained in this case using (135|) 
and (136|1 as, 

g = (37 - mp2 + 2{np2 - l)^^^-^. (39) 

Here it is noticed that when = 0, entry in to the accelerated branch 
is determined by the value of 7. And, 7 = 1/3 forms the critical value that 
separates the accelerating phase from the deccelerating, compared to FRW 
case where 7 = 2/3 is the critical value. 



(a) 




(b)- 



Figure 5. Plot for decceleration parameter q for e = +1 (a) shows different values 
of shear parameter S and Vc = 0.1, (b) shows evolution with different values of Tc 
and S = 0.001 with Q.p^ = 0.04 ,7 = 4/3, K = Q . 
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Figure 6. Figure (a) shows behaviour of $7^2 in e = — 1 case for different values of 
ripg and 7 = 4/3. (b) for different values of 7 and ilp^ = 0.04. In both the plots 
we use S = 0.001, K = and = 0.1. 

4 Conclusions 

We studied the evolution of cosmological density parameter in an anisotropic 
DGP braneworld model. The two branches of solutions are considered sepa- 
rately but we emphasized on the self accelerated branch with shear parameter 
and crossover scale. It is found that the evolution of dynamical density pa- 
rameter is not altered much during the late universe by the presence of shear. 
It is seen in figure 1 that asymptotic value of Qp as S ^ approaches 0.5 
and decreases for large values of S . The behaviour of density parameter 
for different values of shear is shown in figure 2 and is compared with the 
isotropic case. In the early universe shear plays a significant role and in the 
late universe evolution coincides with that of a isotropic DGP case. It is 
also observed that in the presence of shear irrespective of the value 7, the 
behaviour of Qp is same near S = 0. The density parameter Qp does not 
seem to vary much with value of cross over scale, for a smaller value of shear, 
as shown in the figure 4(a). Increasing the value of shear to a higher value 
(S = 0.01) the density parameter varies significantly, depending on cross 
over scale, in the early universe. But in the late universe is insensitive to the 
value of Tc- Deceleration parameter is plotted in figure 5, it shows that in 
the presence of shear universe enters a self accelerating phase at later time, 
for high values of shear. It is observed that value of Vc does not change the 
evolution of q significantly and universe enters a self-accelerated phase as 
shown in figure 5(b). Now, in the e = — 1 or fully five dimensional branch of 
the solution, the density parameter shows interesting behaviour depending 
on the value of 7, in the limit S* — ?■ the evolution is drastically changed for 
7 = 2/3. Looking at the decceleration parameter in this case, we found that 
7 = 1/3 forms the critical value that divides the accelerating phase from 
the deccelerating. In this work we were interested in the K = case for 
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both the branches of DGP model, the behaviour of ilp and q would change 
significantly for the universe with a non-zero spatial curvature. 

As the anisotropic effects are relevant at high redshifts and decreases with 
the expansion of the universe, we would like to estimate the value of the shear 
at the last scattering surface. We obtain the value of shear from our equations 
to be of the order ~ 1.68 x 10^^*^, this is in agreement with the limits 
obtained for the other Bianchi models. Hence, the estimated value of shear 
from our model is compatible with the late time evolution of DGP model 
and the crossover scale in this case is consistent with the one obtained by 
others [HI Ho]. Finally, we would hke to conclude that the dynamical evolution 
of density parameter does not affect late time behaviour of the DGP model 
in the anisotropic case and the results are in agreement with isotropic case 
in the late universe. 
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